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Abstract 

We analyze solvability of a special form of distributed order fractional 
differential equations 

2 4>ih)D~'y(t)dj= [ 2 Mi)D''z(t)d 7 , t>o, (i) 

Jo 

within <S+, the space of tempered distributions supported by [0, oo). 
Keywords: Distributed order fractional differential equations; Tempered distri- 
bution; Laplace transform 
AMS Subject 26A33; 46F12 



1 Motivation and Introduction 

We consider a distributed order fractional differential equation ([T]) which arises 
in the theory of constitutive equations for viscoelastic bodies. 4>i, <f>2 are certain 
functions or distributions which characterize a material under consideration and, 
in general, are determined from experiments. Z? 7 , 7 € 1, is the left Riemann- 
Liouville operator of fractional differentiation or integration defined as follows. 

Denote by Lj oc+ (M) the space of locally integrable functions y on R such 
that y(t) = 0, t < 0. Then for y £ Lj oc+ (M) the left fractional integral of order 
7 > is defined by 

FvW^^rjit-ry-iyWdT, t>0. 

Here T is the Euler gamma function, ff 7 = then I°y :— y. It can be shown 
(cf. [32]) that for y £_L} oc (R) lim 7 ^ ^ 7 y(i) = v(t), t £ E almost everywhere. 

Denote by AC k (R + ) the space of functions y such that y has continuous 
derivatives on R + = {x £ R; x > 0} up to the order k — 1 and k-th derivative 
is locally integrable function. We extend such functions to R so that y(i) = 0, 
t < 0. 
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Let y G AC k (M. + ). Riemann-Liouville's fractional derivative of order 7 > 0, 
7 < k for some k € N, is defined by 

Wy(t) := |^ 7 2/(*), t > 0. 

It follows that D~'y G L} oc+ {R). We refer to Section [2] for the definition of D~'y, 
y G S' + and 7 € E. 

Note that D 1 l 1 y = y for y G L, 1 oc+ (R) and PD^y = y, 7 > in the sense 
of tempered distribution. We sometimes denote D~ 1 y = / 7 y, 7 > 0. 

Let (/> be continuous function in [c,d] C [0, fc], A; G N. Distributed order 
fractional derivative of y G AC k (R+) is given by 

d 

<f>( 7 )Dty(t)d 7 . 

Equation JTJ) models various physical processes. For example, if it models a 
viscoelastic body then fl} represents a constitutive equation of a material and 
connects strain y(t) with corresponding stress z(t) at time instant t > 0. For a 
standard linear viscoelastic body the constitutive equation is given as 

y{t)+byV{t)=z{t)+azV{t), 

where (-)^ = JjO and a, 6 are experimentally determined constants with the 
restriction < a < b following from the second law of thermodynamics. A 
slight generalization (see [2] and references there) of this equation is achieved 
by replacing the first derivative by a derivative of real order a > 

y(t) + bD a y(t) =z(t)+aD a z(t), (2) 

where, again 0<a<6. If < a < 1 then ([2]) represents visco-elastic effects 
while for 1 < a < 2, describes visco- inert ial effects of a material. Standard 
procedure in building rheological models is to use more then one derivative on 
each side of constitutive equation. When this is done (J2J) becomes 

N M 

b n D^y(t) = 2 a m D a ™z(t), (3) 

n— rn—0 

where M, N G N, a m , b„ £ R and a m , (3 n G R, < a m ,(3 n < 2. 

Equation ^ is interpreted in [T] as a Riemann sum. Moreover, in [T] is pro- 
posed the constitutive equation of a linear visco-elastic body in a " distributed" 
order form as 

2 rl 

Ml)D"<y(t)dj = / 2 (7)^z(i)d7, t > 0. (4) 



In model ([U all derivatives of the stress D~> ' z depend on all derivatives of the 
strain D 1 y for 7 G [c, d]. Since the upper bound in integrals in ([U is two, both, 
visco-elastic and visco-inertial effects are included. The presence of integral on 
the left hand side indicates, as experiments show, that dissipation properties 
depend on the order of the derivative. The integral on the right hand side is a 
consequence of the known principle of equipresence. 
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In this paper we are looking for an S' + solution z to ^ for a given but 
arbitrary y G S', . Such solution will be used in [5] for solving a differential 
equation of motion coupled with constitutive equation ((H). 

In Section [2] we extend results obtained in [5] concerning integral in (f¥|). 
Afterwards we define distributed order fractional derivative in S' + and derive 
it's main properties. In Section[3]we state without proof (which is given in [5]) 
a theorem on the existence and uniqueness of a solution to a linear fractional 
differential equation in the frame of S',. Also we derive properties of such 
solution. In SectionU]we connect condition for the uniqueness with a dissipation 
inequality that guarantees physical admissability of a equation in |4]). 

We note that the models with distributed order derivatives were analyzed 
for example in [1 E 13 E3 EE] • 



2 Distributed order fractional derivative 

We denote by 5(R) the space of rapidly decreasing functions in R and by £'(11$) 
its dual, the space of tempered distributions; <S+(R) denotes its subspace con- 
sisting of distributions supported by [0, oo). In the sequel we drop R in the 
notation. We consider in S' + the family 

/.<*> = I ieB '. a>0 ' (•) 



d 

dx 1 



f a+N (x), a < 0, a + N > 0, TV G N, 



where H is Heaviside's function. It is known that f a * fp — fa+p, a, /3 G R. The 
convolution operator f a * in S' + is the operator of fractional differentiation for 
a < and of fractional integration for a > 0. It coincides with the operator of 
derivation for —a G N and integration for a G N. Let a > and y G L 1 1 oc+ (R). 
Then I a y = f a *y. Let y G AC k {W + ) and < a < k. Then D a y = f_ a * y. 
Recall, if y G S' + then its Laplace transform is defined by 

y(s) = Cy(s) = (y(t), ^(t)e- st ), Re s > 0, a G R, 

where <p G C°°, y> = I on (—a, oo) and <p — in (— oo, —2a), a > 0. Note that 
Cy is an analytic function for Re s > and that the definition of Cy does not 
depend on a chosen function ip with given proprieties. We will often use the 
identity 

C{f a *y){s) = -ly(s), Re S >o. 

First we analyze integral J 4>{^)D 1 y{-)d'^. To do this we examine the 

supp tf> 

mapping a i— > : R — > <Sj_, for given y G 5^ (in [S] we have considered 

y eLl c (R)nS' + ). 

Proposition 2.1. (a) Let a G R be fixed. Then the mapping y i— > -D"j/ is linear 
and continuous from S' + to S', . 

(b) Let y G S' + be fixed. Then a t— > Z? Q y is a smooth mapping from R to <S^_. 
(cj T/ie mapping (a,y) -D a y is continuous from R x 5^ to <S^_. 

Proof, (a) The continuity of y i— > Z?"y = /_« * y is clear since for g G 
/ l— * / * 9 is a continuous mapping of into S' + . 
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(b) It is known that there exists a continuous function F, supp F C [0, oo) and 
fc G N such that \F(x)\ < C(l + \x\) k , i£l and y — D k F. So the mapping 
a i— > -D"y equals a i— > D a+k F. By [SI Proposition 1] we know that for fixed fc 
and a € R, a + fc D a+k F is smooth so the same hold for a i— > D a+k F. 

(c) Since 5 is Frechet space as well as locally convex, the separate continuity 
proved in (a) and (b) imply joint continuity (c.f. [131 Corollary to Theorem 
34.1]) . □ 

By £'(R) is denoted the space of compactly supported distributions i.e. the 
dual space of £(R) = Qf (R). 

Definition 2.2. Let (j) G £'(R) and y £ S' + . Then J <f>{i)D^y drf is defined 

SUpp 

as an element of S!±_ by 



( J (t>{ 1 )Diy{t)d 1 ,y{t)) = (ct>{ 1 ),{D 1 y{t),y{t))), V G 5(R). (6) 

supp </> 

5mc/i defined distribution is called distributed order fractional derivative. 

By Proposition 12. 1[ part (b), 7 ► Z3 7 y : R — » Si is smooth as well as 
7 1— > (D 7 y(t),ip(t)) : R — > R. Since 5 is a Frechet space it follows that in 
it's dual space the strong and weak boundedness are the same, thus a linear 
functional defined by © is continuous from S to C and therefore is a tempered 
distribution supported by [0, 00). 

The following two examples are often used in applications. 

Example 2.3. a) Let 7i e R, i G {0,l,..,fc} and <£(•) = £i=o a ^ W G - 7i). 
Then (JSJ) gives 

/ 0( 7 )I7M-)d7 = ^ ai DM-), ^ 5V 

SUpp 

b) Let ^ be a continuous function in [c, d] C R for some c < d, then 

/ «/)(7)^y(-)d7= hm V0(7i)£» 7i y(-)A7i, in 5^, 

where 7j are points of interval [c, d] in usual definition of the Riemann sum for 
the integral. 

Proposition 2.4. Let 4> G £'(R) and y ES' + . Then: 
a) 

J (t>{l)D~<yd"i 

SUpp (j) 

is a linear and continuous mapping from to S'^_ . 
b) 

£( J 0( 7 )^d 7 )(s) =y(a)(0(7).* 7 ). i?e S >0. 

SUpp 
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c) If 4> is continuous function on [c, d] and <p(^f) — for 7 ^ [c, d] then 

rd rd 

£{ / ^)Lnydi){a) = y(s) / 0( 7 )s 7 d 7 , i?es > 0. 



Proof, a) Clearly, this mapping is linear. Let 2M — > in <S^_. Then (0(7), {D 7 y ni ip)} 
0, as n — ► 00, since by Proposition 12.11 part (a), (Z? 7 y„, tp) — > 0, as n — > 00. 

b) By the definition, 

£( / cf>( 7 )D"<yd 1 )(s) = ( J <t>( 7 )D"<y(t)d 7 ,<p(t)e- st ) 

supp 4> supp <fi 

= (0( 7 ), (D-'y (t), p(t)e- st )) = (0( 7 ), s^y(s)}, Re a > 0. 

c) In the case that (A is continuous we have 



(</)(7),s 7 2/(s)} = / </)(7)s 7 y(s)d7 

supp 

and therefore the assertion follows. □ 

If we assume that y 7 z £ S' + in (JH), put = (j> 2 and 5 = f anpp< j >1 4>i(j)D~" ' z 
then the solvability of (J4j> with respect to z reduces to the solvability of 

J ^)D^z = g, geS' + . (7) 

J supp <fc 

3 Linear fractional differential equation in S' + 

Assuming that g G S' + and that <p is of the form as in Example 12.31 equation 
([7|) becomes 

k 

^a. ( L) 7i z = g, in S' + . (8) 

i=0 

We suppose that 7 j 6 [0, 2) such that 70 > 7i > 7 i+i > 7fc, i E {1, ...A: — 1}. 
Theorem 3.1. Equation has a unique solution z S (Si. i/ and only if 

k 

{A ) 0, s S C+ = {s e C ; .Re s > 0}. 

i=0 

The proof is given in [6]. The solution to ||5J) that is obtain in Theorem 13. II 
is given by z = I * g, where 

l(t)=£-\— F ^ )(*), t>0, (9) 

Ei= a * s7, 

is a fundamental solution to ([5]) i.e. solution to yV_ Q a^Z? 7 *?/ = <5. 
The following lemma gives main properties of Z defined by 
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Lemma 3.2. Assume {Aq). Let 7.; £ [0,2) and 70 > 7; > "fi+i > "fk, for all 
i £ {1, — 1}. Let I be defined by |Tpp and l(t) = 0, t < 0. Then: 

(i) I is a locally integrable function in R. 

(wj Moreover, I is absolutely continuous in R, 2/70 — 7fc > 1- 

Proof, (i) Let 7^ = and a k ^ 0. Consider the integral 

t > 0, 



where T = U<=i an d f° r arbitrarily chosen i?>0,0<e<i? and xq > 0, L; 
are given by 

To : {z; Re z — xq; < argz < <po — arcsin -77}; 

R 

Fn z = Re itp , -0o < <t>o < 95 < tt; r 2 : z = fle^, -tt < </j < -0 O < 0; 

r 3 : z = ee tv , —ir < (p < n; T 4 : z = le 11 ; r 5 : z = ie" 17r , x e (e, i?). 

By the Cauchy residue theorem, letting e — > and i? — ► 00, one obtains 

n s j 
i(t) = £ Res{ — ^ -} + lo(t), t>0, (10) 

s=s m ,m=l Z^i=O a i s1 ' 

where s m are poles of the function s 1— > _ fc e — and 



io(*) = - [ 
1" Jo 



00 

e~ st r(s)ds, t > 0, (11) 



with 

Etc «^ 271 + 2 £Ji=<M*y W 7i+7j cos( 7 « - 7 j)tt 



= ^ 2E r0 a^ 7i S in(7^) g > Q _ 



We refer to [1] for similar calculations. 
Let < a < b. Then 

poo 1 

Z Q (t)cft = / (e- SQ -e~ s6 )-r(s)ds. (12) 

Since 7^ = 0, this integral is finite in a neighborhood of s = 0. In a neighborhood 
of s = co we have ~ s70 1 +1 . Thus 70 + 1 > 1 implies that the integral in (TT2"]) 
is finite. Therefore, is locally integrable. By (fT0|) and the fact that Re s m < 
(by (Aq)), we obtain that ^ is locally integrable. 
Let jk > 0. Then 

z = £- 1 (^-)*r- 1 (-^ ) = U*h, (13) 

where 

h^C^i—j-^ -), fj l=ll - lkl i 6 {0,1,..,*}, Res>0, 
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and f lk is defined by |(SJ). Note that (3k = 0. By the first part of the proof, l\ is a 
locally integrable function. Since / 7fe is locally integrable, I is locally integrable 
as the convolution of two locally integrable functions. 

(ii) Let 7^. = 0. Then (TTT|) is finite in a neighborhood of s = 0. In a 
neighborhood of s — oo we have that r(s) ~ -^j-. Since 70 > 1, the integrand in 
(fTTj) is integrable for all i > and (dJJ is finite. Let t > 0. Since \e~ st r(s)\ < 
e~ st °\r(s)\ := g(s) for all t > to and g G i 1 ((0, 00)), by the classical theory 
we obtain that (fTTj) defines a continuous function for t > to. It follows that 
lo and Z (by (fTUf ) are continuous for t > 0. Further on, since <9t(e~ st r(s)) < 
e~ sf °\sr(s)\ :— gi(s), for all t > to and <?x S L 1 (0, 00), we obtain that Iq is 
differentiable and 

'oW = / (-s)e" s *r(s), i > 0. 
Jo 

Since —sr(s) ~ in a neighborhood of s — 00, as in (i) we show that l' Q is a 

locally integrable function. Therefore, the derivative 1' exists and it is a locally 
integrable function. It means that / is absolutely continuous. 

For 7^ > we proceed as in (i) and obtain I — f lk * l\ with l± absolutely 
continuous. Therefore, / is also absolutely continuous. □ 

Remark 3.3. If ji £ [0, 00) and if 70 > p, p G N, then I is continuous in K 
as well as its derivatives up to order p — 1 while the p-th derivative is a locally 
integrable function, i.e. I G AC P . 



4 Comments from mechanics and further appli- 
cations 

Equation ([8]) represents a constitutive equation of a visco-elastic body. We will 
show that in the case when there exist so G C+ such that J"V._ c^Sp 1 = 
it follows that the dissipation inequality (fTTj) . (see 0) is violated. The dis- 
sipation inequality requires that for any T > 0, any y and z the solution to 
Si=o 0'iD li z(t) = y(t), t > 0, the dissipation work, Ad is nonnegative, i.e. 

A d = [ z(t)yW(t)dt>0. (14) 



Let T > and y(t) = H(t) - H(t - r), < r < T, t > 0. Then 

z(t) = C-\-— r ^ 6 ^ ) = f g(u)du, t>0, (15) 

where g(u) = l(u) — l(u — r) and Z is the fundamental solution to ((8j) given by 
©. Assume that u>(sq) = XV=o a i s Q i = ^ f° r s o = u + u > 0. Then by (fTU|) 



wy L )(s) •f— J w { . L '{s) 1 



Further, note that 



e 



[cos(vt) + i sin(i>i)] 



W 



£ > 0, (16) 
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represents oscillations with increasing amplitudes. Inserting (|15p in (114[) we 
obtain 

A d = lim / g(u)du - / g(u)du > 0, T,t > 0. (17) 



t->o 



It is obvious that due to the presence of the term (fT6| in g(i) the inequality (fT7|) 
could be violated by a suitable choice of t. 
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